Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions.

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

Plant MicroRNA Potential in Targeting Sars-CoV-2 Genome Offering Efficient Antiviral MiRNA-Based Therapies.

BACKGROUND: In 2019, severe acute respiratory coronavirus II (or SARS-COV-2) emerged in Wuhan, China, rapidly becoming a global pandemic. Coronavirus genus (Coronaviridae) has the largest single-stranded positive-sense RNA genome (~30 kb) among the human infected single-stranded RNA viruses. OBJECTIVES: For the study of active therapeutic plant-derived miRNA(s), it may be possible to uptake the miRNAs and their biological role in the host cell. In this study, we bioinformatically searched plant miRNAs that can potentially interact with the Sars-CoV-2 genome within the 3'- UTR region and have prompt antiviral activity. MATERIALS AND METHODS: We searched the plant miRNAs that target the 3'-UTR flanking region of the Sars-CoV-2 genome by employing the RNAHybrid, RNA22, and STarMir miRNA/target prediction tools. RESULTS: The RNAHybrid algorithm found 63 plant miRNAs having hybridization energy with less or equal to -25 kcal.mol-1. Besides, RNA22 and STarMir tools identified eight interactions between the plant miRNAs and the targeted RNA sequence. pvu-miR159a. 2 and sbi-miR5387b were predicted as the most effectively interacting miRNAs in targeting the 3'-UTR sequence, not only by the RNA22 tool but also by the STarMir tool at the same position. However, the GC content of the pvumiR159a. 2 is 55% instead of sbi-miR5387b, which is a GC enriched sequence (71.43%) that may activate TLR receptors. CONCLUSION: In our opinion, they are potent plant-derived miRNA candidates that have a great chance of targeting the Sars-CoV-2 genome in the 3'-UTR region in vitro. Therefore, we propose pvu-miR159a.2 for studying antiviral miRNA-based therapies without any essential side effects in vivo.

Highly conserved s2m element of SARS-CoV-2 dimerizes via a kissing complex and interacts with host miRNA-1307-3p.

The ongoing COVID-19 pandemic highlights the necessity for a more fundamental understanding of the coronavirus life cycle. The causative agent of the disease, SARS-CoV-2, is being studied extensively from a structural standpoint in order to gain insight into key molecular mechanisms required for its survival. Contained within the untranslated regions of the SARS-CoV-2 genome are various conserved stem-loop elements that are believed to function in RNA replication, viral protein translation, and discontinuous transcription. While the majority of these regions are variable in sequence, a 41-nucleotide s2m element within the genome 3' untranslated region is highly conserved among coronaviruses and three other viral families. In this study, we demonstrate that the SARS-CoV-2 s2m element dimerizes by forming an intermediate homodimeric kissing complex structure that is subsequently converted to a thermodynamically stable duplex conformation. This process is aided by the viral nucleocapsid protein, potentially indicating a role in mediating genome dimerization. Furthermore, we demonstrate that the s2m element interacts with multiple copies of host cellular microRNA (miRNA) 1307-3p. Taken together, our results highlight the potential significance of the dimer structures formed by the s2m element in key biological processes and implicate the motif as a possible therapeutic drug target for COVID-19 and other coronavirus-related diseases.

Deciphering the link between Diabetes mellitus and SARS-CoV-2 infection through differential targeting of microRNAs in the human pancreas.

PURPOSE: Coronavirus Disease 2019 (COVID-19) severity and Diabetes mellitus affect each other bidirectionally. However, the cause of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection on the incidence of diabetes is unclear. In the SARS-CoV-2-infected cells, host microRNAs (miRNAs) may target the native gene transcripts as well as the viral genomic and subgenomic RNAs. Here, we investigated the role of miRNAs in linking Diabetes to SARS-CoV-2 infection in the human pancreas. METHODS: Differential gene expression and disease enrichment analyses were performed on an RNA-Seq dataset of human embryonic stem cell-derived (hESC) mock-infected and SARS-CoV-2-infected pancreatic organoids to obtain the dysregulated Diabetes-associated genes. The miRNA target prediction for the Diabetes-associated gene transcripts and the SARS-CoV-2 RNAs has been made to determine the common miRNAs targeting them. Minimum Free Energy (MFE) analysis was done to identify the miRNAs, preferably targeting SARS-CoV-2 RNAs over the Diabetes-associated gene transcripts. RESULTS: The gene expression and disease enrichment analyses of the RNA-Seq data have revealed five biomarker genes, i.e., CP, SOCS3, AGT, PSMB8 and CFB that are associated with Diabetes and get significantly upregulated in the pancreas following SARS-CoV-2-infection. Four miRNAs, i.e., hsa-miR-298, hsa-miR-3925-5p, hsa-miR-4691-3p and hsa-miR-5196-5p, showed preferential targeting of the SARS-CoV-2 genome over the cell's Diabetes-associated messenger RNAs (mRNAs) in the human pancreas. CONCLUSION: Our study proposes that the differential targeting of the Diabetes-associated host genes by the miRNAs may lead to diabetic complications or new-onset Diabetes that can worsen the condition of COVID-19 patients.

Potential Therapeutic Effect of Micrornas in Extracellular Vesicles from Mesenchymal Stem Cells against SARS-CoV-2.

Extracellular vesicles (EVs) are cell-released, nanometer-scaled, membrane-bound materials and contain diverse contents including proteins, small peptides, and nucleic acids. Once released, EVs can alter the microenvironment and regulate a myriad of cellular physiology components, including cell-cell communication, proliferation, differentiation, and immune responses against viral infection. Among the cargoes in the vesicles, small non-coding micro-RNAs (miRNAs) have received attention in that they can regulate the expression of a variety of human genes as well as external viral genes via binding to the complementary mRNAs. In this study, we tested the potential of EVs as therapeutic agents for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. First, we found that the mesenchymal stem-cell-derived EVs (MSC-EVs) enabled the rescue of the cytopathic effect of SARS-CoV-2 virus and the suppression of proinflammatory responses in the infected cells by inhibiting the viral replication. We found that these anti-viral responses were mediated by 17 miRNAs matching the rarely mutated, conserved 3'-untranslated regions (UTR) of the viral genome. The top five miRNAs highly expressed in the MSC-EVs, miR-92a-3p, miR-26a-5p, miR-23a-3p, miR-103a-3p, and miR-181a-5p, were tested. They were bound to the complemented sequence which led to the recovery of the cytopathic effects. These findings suggest that the MSC-EVs are a potential candidate for multiple variants of anti-SARS-CoV-2.

Identifying potential human and medicinal plant microRNAs against SARS-CoV-2 3'UTR region: A computational genomics assessment.

The coronavirus disease of 2019 (COVID-19) began as an outbreak and has taken a toll on human lives. The current pandemic requires scientific attention; hence we designed a systematic computational workflow to identify the cellular microRNAs (miRNAs) from human host possessing the capability to target and silence 3'UTR of SARS-CoV-2 genome. Based on this viewpoint, we extended our miRNA search to medicinal plants like Ocimum tenuiflorum, Zingiber officinale and Piper nigrum, which are well-known to possess antiviral properties, and are often consumed raw or as herbal decoctions. Such an approach, that makes use of miRNA of one species to interact and silence genes of another species including viruses is broadly categorized as cross-kingdom interactions. As a part of our genomics study on host-virus-plant interaction, we identified one unique 3'UTR conserved site 'GGAAGAG' amongst 5024 globally submitted SARS-CoV-2 complete genomes, which can be targeted by the human miRNA 'hsa-miR-1236-3p' and by Z. officinale miRNA 'zof-miR2673b'. Additionally, we also predicted that the members of miR477 family commonly found in these three plant genomes possess an inherent potential to silence viral genome RNA and facilitate antiviral defense against SARS-CoV-2 infection. In conclusion, this study reveals a universal site in the SARS-CoV-2 genome that may be crucial for targeted therapeutics to cure COVID-19.

SNPs at 3'UTR of APOL1 and miR-6741-3p target sites associated with kidney diseases more susceptible to SARS-COV-2 infection: in silco and in vitro studies.

Acute Kidney Injury (AKI) is a common manifestation of COVID-19 and several cases have been reported in the setting of the high-risk APOL1 genotype (common genetic variants). This increases the likelihood that African American people with the high-risk genotype APOL1 are at increased risk for kidney disease in the COVID-19 environment. Single-nucleotide polymorphisms (SNPs) are found in various microRNAs (miRNAs) and target genes change the miRNA activity that leads to different diseases. Evidence has shown that SNPs increase/decrease the effectiveness of the interaction between miRNAs and disease-related target genes. The aim of this study is not only to identify miRSNPs on the APOL1 gene and SNPs in miRNA genes targeting 3'UTR but also to evaluate the effect of these gene variations in kidney patients and their association with SARS-COV-2 infection. In 3'UTR of the APOL1 gene, we detected 96 miRNA binding sites and 35 different SNPs with 10 different online software in the binding sites of the miRNA (in silico). Also we studied gene expression of patients and control samples by using qRT-PCR (in vitro). In silico study, the binding site of miR-6741-3p on APOL1 has two SNPs (rs1288875001, G > C; rs1452517383, A > C) on APOL1 3'UTR, and its genomic sequence is the same nucleotide as rs1288875001. Similarly, two other SNPs (rs1142591, T > A; rs376326225, G > A) were identified in the binding sites of miR-6741-3p at the first position. Here, the miRSNP (rs1288875001) in APOL1 3'UTR and SNP (rs376326225) in the miR-6741-3p genomic sequence are cross-matched in the same binding region. In vitro study, the relative expression levels were calculated by the 2-ΔΔCt method & Mann-Whitney U test. The expression of APOL1 gene was different in chronic kidney patients along with COVID-19. By these results, APOL1 expression was found lower in patients than healthy (p < 0.05) in kidney patients along with COVID-19. In addition, miR-6741-3p targets many APOL1-related genes (TLR7, SLC6A19, IL-6,10,18, chemokine (C-C motif) ligand 5, SWT1, NFYB, BRF1, HES2, NFYB, MED12L, MAFG, GTF2H5, TRAF3, angiotensin II receptor-associated protein, PRSS23) by evaluating online software in the binding sites of the miR-6741-3p. miR-6741-3p has not previously shown any association with kidney diseases and SARS-COV-2 infection. It assures that APOL1 can have a significant consequence in kidney-associated diseases by different pathways. Henceforth, this study represents and demonstrates an effective association between miR-6741-3p and kidney diseases, i.e., collapsing glomerulopathy, chronic kidney disease (CKD), acute kidney injury (AKI), and tubulointerstitial lesions susceptibility to SARS-COV-2 infection via in silico and in vitro exploration and recommended to have better insight.

De novo 3D models of SARS-CoV-2 RNA elements from consensus experimental secondary structures.

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta's FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5' UTR; the reverse complement of the 5' UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3' UTR. For eleven of these elements (the stems in SL1-8, reverse complement of SL1-4, FSE, s2m and 3' UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets ('FARFAR2-SARS-CoV-2', https://github.com/DasLab/FARFAR2-SARS-CoV-2; and 'FARFAR2-Apo-Riboswitch', at https://github.com/DasLab/FARFAR2-Apo-Riboswitch') include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules.

Genome-scale deconvolution of RNA structure ensembles.

RNA structure heterogeneity is a major challenge when querying RNA structures with chemical probing. We introduce DRACO, an algorithm for the deconvolution of coexisting RNA conformations from mutational profiling experiments. Analysis of the SARS-CoV-2 genome using dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) and DRACO, identifies multiple regions that fold into two mutually exclusive conformations, including a conserved structural switch in the 3' untranslated region. This work may open the way to dissecting the heterogeneity of the RNA structurome.

Common Genetic Variation in Humans Impacts In Vitro Susceptibility to SARS-CoV-2 Infection.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant interindividual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly afflict healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants influence vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single-nucleotide polymorphism (rs4702), common in the population and located in the 3' UTR of the protease FURIN, influences alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can have an impact on viral infection and thus contribute to clinical heterogeneity in COVID-19. Ongoing genetic studies will help to identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs.

Tackling the COVID-19 "cytokine storm" with microRNA mimics directly targeting the 3'UTR of pro-inflammatory mRNAs.

COVID-19 is characterized by two major clinical phases, the SARS-CoV-2 infection of target cells and tissues, and a deep inflammatory state, known as "cytokine storm", caused by activation of pro-inflammatory genes, such as NF-kB, STAT-3, IL-6, IL-8, IL-1ß. Among possible anti-inflammatory agents, the "microRNA targeting" should be carefully considered, since it is well known that microRNAs are deeply involved in the expression of cytokines, chemokines and growth factors. The working general hypothesis is that targeting the microRNA network might be important for the development of therapeutic approaches to counteract the COVID-19 induction of inflammatory response. This hypothesis is based on several publications demonstrating the use of miRNA mimics for inhibitory effects on the production of proteins characterizing the COVID-19 "cytokine storm".

The RNA Architecture of the SARS-CoV-2 3'-Untranslated Region.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic. The 3' untranslated region (UTR) of this ß-CoV contains essential cis-acting RNA elements for the viral genome transcription and replication. These elements include an equilibrium between an extended bulged stem-loop (BSL) and a pseudoknot. The existence of such an equilibrium is supported by reverse genetic studies and phylogenetic covariation analysis and is further proposed as a molecular switch essential for the control of the viral RNA polymerase binding. Here, we report the SARS-CoV-2 3' UTR structures in cells that transcribe the viral UTRs harbored in a minigene plasmid and isolated infectious virions using a chemical probing technique, namely dimethyl sulfate (DMS)-mutational profiling with sequencing (MaPseq). Interestingly, the putative pseudoknotted conformation was not observed, indicating that its abundance in our systems is low in the absence of the viral nonstructural proteins (nsps). Similarly, our results also suggest that another functional cis-acting element, the three-helix junction, cannot stably form. The overall architectures of the viral 3' UTRs in the infectious virions and the minigene-transfected cells are almost identical.

Secondary structure determination of conserved SARS-CoV-2 RNA elements by NMR spectroscopy.

The current pandemic situation caused by the Betacoronavirus SARS-CoV-2 (SCoV2) highlights the need for coordinated research to combat COVID-19. A particularly important aspect is the development of medication. In addition to viral proteins, structured RNA elements represent a potent alternative as drug targets. The search for drugs that target RNA requires their high-resolution structural characterization. Using nuclear magnetic resonance (NMR) spectroscopy, a worldwide consortium of NMR researchers aims to characterize potential RNA drug targets of SCoV2. Here, we report the characterization of 15 conserved RNA elements located at the 5' end, the ribosomal frameshift segment and the 3'-untranslated region (3'-UTR) of the SCoV2 genome, their large-scale production and NMR-based secondary structure determination. The NMR data are corroborated with secondary structure probing by DMS footprinting experiments. The close agreement of NMR secondary structure determination of isolated RNA elements with DMS footprinting and NMR performed on larger RNA regions shows that the secondary structure elements fold independently. The NMR data reported here provide the basis for NMR investigations of RNA function, RNA interactions with viral and host proteins and screening campaigns to identify potential RNA binders for pharmaceutical intervention.

Global cataloguing of variations in untranslated regions of viral genome and prediction of key host RNA binding protein-microRNA interactions modulating genome stability in SARS-CoV-2.

BACKGROUND: The world is going through the critical phase of COVID-19 pandemic, caused by human coronavirus, SARS-CoV-2. Worldwide concerted effort to identify viral genomic changes across different sub-types has identified several strong changes in the coding region. However, there have not been many studies focusing on the variations in the 5' and 3' untranslated regions and their consequences. Considering the possible importance of these regions in host mediated regulation of viral RNA genome, we wanted to explore the phenomenon. METHODS: To have an idea of the global changes in 5' and 3'-UTR sequences, we downloaded 8595 complete and high-coverage SARS-CoV-2 genome sequence information from human host in FASTA format from Global Initiative on Sharing All Influenza Data (GISAID) from 15 different geographical regions. Next, we aligned them using Clustal Omega software and investigated the UTR variants. We also looked at the putative host RNA binding protein (RBP) and microRNA binding sites in these regions by 'RBPmap' and 'RNA22 v2' respectively. Expression status of selected RBPs and microRNAs were checked in lungs tissue. RESULTS: We identified 28 unique variants in SARS-CoV-2 UTR region based on a minimum variant percentage cut-off of 0.5. Along with 241C>T change the important 5'-UTR change identified was 187A>G, while 29734G>C, 29742G>A/T and 29774C>T were the most familiar variants of 3'UTR among most of the continents. Furthermore, we found that despite the variations in the UTR regions, binding of host RBP to them remains mostly unaltered, which further influenced the functioning of specific miRNAs. CONCLUSION: Our results, shows for the first time in SARS-Cov-2 infection, a possible cross-talk between host RBPs-miRNAs and viral UTR variants, which ultimately could explain the mechanism of escaping host RNA decay machinery by the virus. The knowledge might be helpful in developing anti-viral compounds in future.

miR-142a-3p promotes the proliferation of porcine hemagglutinating encephalomyelitis virus by targeting Rab3a.

Porcine hemagglutinating encephalomyelitis virus (PHEV) is a typical neurotropic coronavirus that mainly invades the central nervous system (CNS) in piglets and causes vomiting and wasting disease. Emerging evidence suggests that PHEV alters microRNA (miRNA) expression profiles, and miRNA has also been postulated to be involved in its pathogenesis, but the mechanisms underlying this process have not been fully explored. In this study, we found that PHEV infection upregulates miR-142a-3p RNA expression in N2a cells and in the CNS of mice. Downregulation of miR-142a-3p by an miRNA inhibitor led to a significant repression of viral proliferation, implying that it acts as a positive regulator of PHEV proliferation. Using a dual-luciferase reporter assay, miR-142a-3p was found to bind directly bound to the 3' untranslated region (3'UTR) of Rab3a mRNA and downregulate its expression. Knockdown of Rab3a expression by transfection with an miR-142a-3p mimic or Rab3a siRNA significantly increased PHEV replication in N2a cells. Conversely, the use of an miR-142a-3p inhibitor or overexpression of Rab3a resulted in a marked restriction of viral production at both the mRNA and protein level. Our data demonstrate that miR-142a-3p promotes PHEV proliferation by directly targeting Rab3a mRNA, and this provides new insights into the mechanisms of PHEV-related pathogenesis and virus-host interactions.